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Specific heat of the layered Na^CoC^ (x=0.65, 0.70 and 0.75) oxides has been measured in the 
temperature range of 3-360 K and magnetic field of and 9 T. The analysis of data, assuming 
the combined effect of inter-layer superexchange and the phase separation into mesoscopic magnetic 
domains with localized spins embedded in a matrix with itinerant electronic character, suggests that 
the dominant contribution to the specific heat in the region of short-range ordering is mediated by 
quasi-2D antiferromagnetic clusters, perpendicular to the C0O2 layers. 



PACS numbers: 71.27.+a, 73.43.Nq, 74.20.Mn 

The recently discovered superconductivity 1 , and the 
unique combination of high thermopower with low 
electrical resistivity 2 brought the sodium cobalt oxide 
Na^CoC^ and its hydrated counterpart into the focus 
of interest of the last years. Theoretical considerations 3 
point at the essential role of magnetic fluctuations in both 
the superconducting pairing mechanism and the elec- 
tronic transport. The character of fluctuations is closely 
associated with the magnetic dimensionality and frus- 
tration, consequently, capturing the nature of magnetic 
interactions in this system in the widest possible context 
appears to be crucial. 

The Naa;Co02 system is considered a quasi-2D doped 
Mott insulator with antiferromagnetic (AF) interaction 
between the Co atoms on triangular lattice. Sodium be- 
tween the close-packed C0O2 planes serve as a charge 
reservoir, by donating one electron converts the origi- 
nally magnetic Co 4+ (3d 5 ) ions with spin S=l/2 into 
nonmagnetic Co 3+ (3d 6 ) ones. The transport and mag- 
netic properties are therefore very sensitive to the dop- 
ing level and to the degree of frustration, introduced by 
triangular geometry. The system can be doped over a 
wide range, and for most x, Na^OoC^ is metallic. How- 
ever in a narrow range around x =0.5, it represents a 
charge ordered AF insulator—. Charge ordering as a pos- 
sible result of frustration and strong correlations has been 
predicted within several theoretical considerations^* 6 - for 
other doping levels as well. Experimental indications 
of this process were observed by nuclear magnetic res- 
onance studies^. Phase separation on mesoscopic scale 
into Co 4+ -rich magnetic domains and Co 3+ -rich regions 
with itinerant electronic character was strongly indicated 
by resonance methods^ and neutron scattering^. Mag- 
netic phase transition has been observed only for the 
range 0.7<x<0.95 with clear AF spin density wave re- 
solved for 3;>0.8i2iii. Theoretical first principle calcula- 
tions based on local density approximations (LDA) are 



not conclusive concerning the nature of ground state 
and type of interactions. More specifically, the LDA 
calculations 1 ^ predict an itinerant ferromagnetic (FM) 
ground state, whereas LDA+U 5 . calculations admit also 
AF ground states. The corresponding calculations for 
tripled supercells for doping £=2/3 lead to solutions fa- 
voring AF stacking of FM planes 14 . Analysis of the ex- 
change paths in the system has revealed the feasibility 
of strong inter-planar superexchange, mediated by by p z 
orbitals of O and considerably facilitated by empty sp 2 
hybrid orbitals of Na. Experimentally, although the sus- 
ceptibility behavior indicated antiferromagnetism, strong 
in-plane FM fluctuations have been identified by inelastic 
neutron scattering^- 5 ". Consistently with the anisotropic 
character of the structure in which the inter-plane Co dis- 
tances are almost 4 times larger than the in-plane ones, 
the transport was found to be highly anisotropic 16 and 
naturally a quasi-2D magnetic structure has been pre- 
sumed. Nevertheless, the scheme concerning the mag- 
netic dimensionality is recently being changed: an in- 
tense debate brisked up about the possible 3D character 
of magnetism in this structurally layered system. Recent 
polarized- and unpolarized-neutron scattering studyS* 1 ^ 
suggested surprisingly almost isotropic 3D magnetic in- 
teractions, specifically, FM exchange within the layers 
with exchange constant ~ 6 meV and AF correlation be- 
tween the C0O2 layers with ~ 12 meV. 

Specific heat on this system has been predominately 
analysed using the Debye model and the physical in- 
terpretations concerned the discussions about the Som- 
merfeld coefficienliiliiiiii. We present an alternative ap- 
proach, based on mesoscopic phase separation, taking 
into account the most recent neutron scattering results. 
The problem of how will the proposed interactions and 
the existence of phase separation manifest themselves in 
in the specific heat, is the main subject of the present 
work. 
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Our Na x Co02 powder samples with nominal Na con- 
tent of £=0.65, 0.70, 0.75 (hereafter denoted as Na65, 
Na70 and Na75, respectively) in the doping regime, 
where magnetic order is about to form, were prepared 
by the rapid heat-up method^. X-ray powder diffrac- 
tion confirmed the samples to be single phase of hexago- 
nal 7-Naa.Co02 with lattice parameter a around 2.827 
A for all samples, while c=10.939, 10.907 and 10.892 
A for x=0.65, 0.70 and 0.75, respectively. Specific heat 
measurements from 3 K to 360 K have been performed 
in zero magnetic field and in magnetic field of 9 T using 
the conventional Quantum Design PPMS-9 equipment. 

Fig.l shows the temperature behavior of specific heat 
for the Na65 sample, this behavior is similar for all sam- 
ples. The lattice contribution has been estimated by 
well established procedure based on harmonic approx- 
imation of the phonon spectrum, using both the De- 
bye and Einstein model and involving corrections for 
anharmonicityi^iSI. In the fitting procedure the simple 
Sommerfeld term C e =j.T has been used to approximate 
the electronic part of specific heat, yielding the value of 
7=29 mJ/molK 2 . This value is in good agreement with 
that reported in the literaturei2*ii*i£. After separation of 
the electronic and phonon background three features be- 
came pronounced in the residual specific heat (Fig.l right 
axis): 1) There is a broad maximum around 250 K. Ac- 
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FIG. 1: Squares: specific heat experimental data (C) for 
the Na65 sample. Dashed line: fit of the electronic (C e ) and 
phonon (C p hon) contribution. Circles: residual specific heat 
after subtracting the electronic and phonon contribution. 

cording to the literature^, this is the temperature range 
where charge ordering occurs, thus it might be speculated 
that this bump is a possible indicium of charge order- 
ing in specific heat. The entropy change accompanying 
this process we found to be proportional to the electron 
doping. Specific heat measurements on single crystals 
also have evidenced two sharp peaks in this temperature 
rcgionii that could be connected to the anomaly observed 
in our measurements, considerably broadened because of 
the powder state of the samples. 

2) The broad maximum around 40-50 K is attributed to 
short-range ordering (SRO) in the magnetic domains, as 



it will be discussed later. 

3) There is a rather smeared phase transition at T* ~ 28 
K with peak position and shape same in magnetic field 
of and 9 T (Fig. 2). The fact that the transition is not 
rounded by magnetic field could imply that the main in- 
teraction in this compound is AF. This is in agreement 
with the Curie- Weiss fits of the susceptibility, that re- 
vealed a negative Weiss temperature with values around 
120 K for all three samples^. The transition temperature 
is independent of Na content x, on the other hand, the 
area of the peak is suppressed by doping. This is consis- 
tent with the phase separation scenario and suggests that 
the phase, the transition is related to, does not change, 
only its volume fraction decreases with doping. 
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FIG. 2: The experimental C/T vs. T data measured in zero 
magnetic field and in the magnetic field of 9 T. The Na70 
and Na75 data are shifted for clarity down by 0.01 and 0.02 
J/molK 2 , respectively. 

Now, let us get closer look at the magnetic part of the 
specific heat of these unconventional oxides, with the aim 
to interpret the maximum ascribed to SRO. In this anal- 
ysis, keeping in mind the suppressing role of Na on the 
number of magnetic ions, within a simple model only 35% 
(30%, 25%) of Co atoms in the 65% (70%, 75%) electron- 
doped system are considered to be magnetic. The vol- 
ume fraction of magnetic ions of about 21%, estimated 
by muon spin rotation in Nao.ysCoO^, shows that the 
above numbers are reasonable. In the subsequent discus- 
sion the specific heat data will be related to one mole 
of magnetic Co 4+ ions. Assuming phase separation, we 
examine the applicability of the predictions of localized 
AF Heisenberg models to the specific heat contribution 
of the magnetic subsystem (i.e. the domains with local- 
ized spins). The choice of Heisenberg model is justified 
by EPR study and by the obtained value of the g-factor 
that generally reflects the anisotropy of magnetic ions. 
More specifically, the estimated value of the g-factor was 
very close to 2 and the asymmetry of the spectra was 
found only slight. It should be stressed that the same 
approach was used in analysis of inelastic neutron scat- 
tering datc&iS. 

The comparison of the rescaled data to different low- 
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dimensional models can be seen in Fig. 3. The first model 
we tested is the triangular lattice model of spins S=l/2 Zi , 
as the magnetic layers were originally expected to cor- 
respond with the structural C0O2 layers. The discrep- 
ancy between the model and data is considerable. For 
all samples the qualitatively best agreement is achieved 
by fitting the quasi-2D square lattice model, involving 
a weak inter-planar coupling J ± — and this is consistent 
with the idea of "magnetic layers" perpendicular to the 
structural ones, as it will be discussed later. However, 
the obtained parameters should be interpreted with cau- 
tion, as the model ascribes the peak to long-range or- 
dering (LRO) induced by interplanar exchange coupling. 
For the obtained value of this parameter the magnetic 
field of 9 T should induce considerable shift of the peak. 
The fact, that no shift is observable, suggests that the 
peak does not reflect LRO mediated by exchange inter- 
actions. The estimated in-plane exchange constant, how- 
ever, is in good agreement with the Weiss temperature 
obtained from the susceptibility. As the precise exchange 
paths along the perpendicular axis are not known and 
on principle the diagonal exchange between the Co ions 
in the neighboring structural layers cannot be excluded, 
we tested also the frustrated square lattice model^, into 
which the next nearest interaction J2 along the diag- 
onal of the plaquette of the lattice is included. The 
broad maximum is quantitatively better described by this 
model, the fitting analysis suggests rather strong frustra- 
tion with J2 I J \= 0.7. For completness, also the AF 
chain models is included, this fit is again rather close to 
the data. 

The above analysis indicates that SRO represents a 
considerable contribution to specific heat and in this way 
the system does not behave like a 3D magnet. The most 
interesting result of the above comparisons is the qual- 
itative agreement with the square lattice models and a 
question arises, whether it is justifiable. Let us put the 
facts together: 

1) The system is separated into mesoscopic domains 
with localized spins and regions with itinerant electronic 
states. 

2) From the aforecited neutron scattering study of Helme 
et al£ follows, that the main AF interaction propagates 
along the c-axis, perpendicular to the C0O2 layers, while 
FM fluctuations are observed within the layers at 6K. It 
must be noted, that in this study the localized 3D spin 
wave approach has been used to analyze the dispersion 
modes, which does not consider the charge ordering or 
phase separation in the system. The observed spin ex- 
citation spectrum however, hints to the possible phase 
separation, namely into FM in-plane clusters of Co 4+ 
ions in a matrix of non-magnetic Co 3+ . Helme et al. 
speculate, that in order to obtain consistency with the 
observed sharp spin modes along the c-axis these clus- 
ters would have to be aligned vertically above each other 
over many layers. 

3) The theory of sodium ordering proposed by Zhang 
et aliSi found stable ordering patterns for different dop- 
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FIG. 3: Comparison of the magnetic part of specific heat 
(symbols) with different low-dimensional AF Heisenberg mod- 
els. 



ings in 7-Na x Co02. Regarding the interplay between Na 
ordering and charge ordering of Co in the planes, consid- 
erably weakened in-plane exchange interactions can be 
expected due to the increased distance between the mag- 
netic Co ions in some specific directions, all the more as 
the magnetic orbitals are oriented out of plane. In other 
words, more or less continuous stripe- like regions of mag- 
netic ions are formed within the layers, depending on Na 
doping level. 

Summing up these facts, the existence of a low dimen- 
sional magnet, most likely in a squared arrangement of 
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magnetic moments with AF interactions perpendicular 
to the C0O2 planes becomes possible. It is notewor- 
thy, that the AF exchange constants obtained from our 
rough fits (120 - 140 K ~ 10 - 12 meV) are in excellent 
agreement with that estimated from neutron scattering 
results^. Moreover, the increasing tendency of this con- 
stant with Na doping supports the reasoning of Johannes 
et ali^ about the superexchange enhanced by Na. It re- 
mains a question, what is the origin of FM correlations. 
One of the possible explanation could be the itinerant 
electrons, nevertheless, this question should be investi- 
gated further. 

In summary, low temperature specific heat analysis 
has shown that the dominant contribution to the specific 
heat in the region of short-range ordering is mediated 
by quasi-2D antiferromagnetic clusters, perpendicular to 
the C0O2 layers in the Na^OoC^ (x=0.65-0.75) system. 
Comparison of the magnetic part of specific heat to the 
predictions of different AF Heisenberg models suggests, 
that the main magnetic interactions are realized on a 
square lattice rather than on a triangular one. The esti- 
mated AF exchange constant is in good agreement with 



the neutron scattering results of Helme et al. It would be 
instructive to study the contribution of the other com- 
ponent of the phase separated system, i.e. the itinerant 
electrons as a possible origin of the observed FM fluctu- 
ations that might become more visible below the phase 
transition. 
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